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Abstract 

The Q-state Potts model can be extended to non-integer and even com- 
plex Q by expressing the partition function in the Fortuin-Kasteleyn (F-K) 
representation. In the F-K representation the partition function, Z(Q,a), 
is a polynomial in Q and v = a — 1 (a = e@ J ) and the coefficients of this 
polynomial, <£(&, c), are the number of graphs on the lattice consisting of b 
bonds and c connected clusters. We introduce the random-cluster transfer 
matrix to compute &(b,c) exactly on finite square lattices with several types 
of boundary conditions. Given the F-K representation of the partition func- 
tion we begin by studying the critical Potts model Zqp = Z(Q, a c (Q)), where 
a c(Q) = 1 + \[Q- We find a set of zeros in the complex w = \[Q plane that 
map to (or close to) the Beraha numbers for real positive Q. We also identify 
Q C (L), the value of Q for a lattice of width L above which the locus of zeros 
in the complex p = v/y/Q plane lies on the unit circle. By finite-size scaling 
we find that 1/Q C (L) — » as L — > oo. We then study zeros of the antiferro- 
magnetic (AF) Potts model in the complex Q plane and determine Q c {a), the 
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largest value of Q for a fixed value of a below which there is AF order. We 
find excellent agreement with Baxter's conjecture Q^ F (a) = (1 — a) (a + 3). 
We also investigate the locus of zeros of the ferromagnetic Potts model in the 
complex Q plane and confirm that Q FM (a) = (a— l) 2 . We show that the edge 
singularity in the complex Q plane approaches Q c as Q C {L) ~ Q c +AL~ y i, and 
determine the scaling exponent y q for several values of Q. Finally, by finite 
size scaling of the Fisher zeros near the antiferromagnetic critical point we 
determine the thermal exponent yt as a function of Q in the range 2 < Q < 3. 
Using data for lattices of size 3 < L < 8 we find that y t is a smooth function 
of Q and is well fit by y t = 1+ ^+f" 2 where u = -f cos" 1 For Q = 3 we 
find y t ~ 0.6; however if we include lattices up to L = 12 we find y t ~ 0.50(8) 
in rough agreement with a recent result of Ferreira and Sokal [J. Stat. Phys. 
96, 461 (1999)]. 
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I. INTRODUCTION 



The Q-state Potts model [|I| in two dimensions exhibits a rich variety of critical behavior 
and is very fertile ground for the analytical and numerical investigation of first- and second- 
order phase transitions. With the exception of the Q = 2 Potts (Ising) model in the absence 
of an external magnetic field, exact solutions for arbitrary Q are not known. However, 
some exact results at the critical temperature have been established for the Q-st&te Potts 
model. From the duality relation the ferromagnetic critical temperature is known to be 
T c = J/k B \n(l + ^fQ) for the isotropic square lattice. Baxter M calculated the free energy 
at T c in the thermodynamic limit, and showed that the Potts model has a second-order 
phase transition for Q < 4 and a first-order transition for Q > 4. The critical exponents for 
the ferromagnetic Potts model are well known 0-0] . 

On the other hand, the antiferromagnetic Potts model is much less well understood than 
the ferromagnetic model. Recently the three-state Potts antiferromagnet on the square 
lattice has attracted a good deal of interest |6|-p6|]. Baxter [|Kj conjectured that the critical 



point of the Potts antiferromagnet on the square lattice is given by T c = J / kB\n(y/A — Q — l), 
and evaluated the critical free energy and internal energy. The Baxter formula for the critical 
temperature gives the known exact value for Q = 2, a critical point at zero temperature 
for Q = 3, and no critical point for Q > 3. For continuous Q in the range < Q < 3, 
Kim et al. |27| have studied the antiferromagnetic Potts critical point through the zeros 
of the partition function and found good agreement with the Baxter formula. With the 
exception of the Ising model, the critical exponents of the Potts antiferromagnets are not 



known. However, for Q = 3 the ratio of critical exponents j/u is known to be 5/3 ||iq , |T7 |. 

By introducing the concept of the zeros of the partition function in the complex magnetic- 
field plane (Yang-Lee zeros), Yang and Lee [^] proposed a mechanism for the occurrence of 
phase transitions in the thermodynamic limit and yielded a new insight into the unsolved 
problem of the Ising model in an arbitrary nonzero external magnetic field. It has been shown 
501 that the distribution of the zeros of a model determines its critical behavior. Lee 
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and Yang |28fl also formulated the celebrated circle theorem which states that the Yang-Lee 
zeros of the Ising ferromagnet lie on the unit circle in the complex magnetic-field (x = e 1311 ) 
plane. However, for the Q-state Potts model with Q > 2 the Yang-Lee zeros lie close to, 
but not on, the unit circle with the two exceptions of the critical point x — 1 (h — 0) itself 



and the zeros in the limit T = |JT| . 



Fisher [|32| emphasized that the partition function zeros in the complex temperature 



plane (Fisher zeros) are also very useful in understanding phase transitions, and showed 
that for the square lattice Ising model in the absence of an external magnetic field the 
Fisher zeros lie on two circles in the thermodynamic limit. In particular, using the Fisher 
zeros both the ferromagnetic phase and the antiferromagnetic phase can be considered at 
the same time. The critical behavior of the Potts model in both the ferromagnetic and 
antiferromagnetic phases have been studied using the distribution of the Fisher zeros, and 



the Baxter conjecture for the antiferromagnetic critical temperature has been verified |27 
Recently the Fisher zeros of the Q-state Potts model on square lattices have been studied 
extensively for integer Q > 2 and noninteger Q |27[| . Exact numerical studies have 



shown P?| , pB] , pDl , ^n| , ^T| , ^5| , P| , ^D| that for self-dual boundary conditions the Fisher zeros of the 
Q > 1 Potts models on a finite square lattice are located on the unit circle in the complex 
p plane for Re(p) > 0, where p = [eP 3 — l)/y/Q. It has been analytically shown that all 
the Fisher zeros of the infinite-state Potts model lie on the unit circle for any size of square 



lattice with self-dual boundary conditions [42], and the Fisher zeros near the ferromagnetic 



critical point of the Q > 4 Potts models on the square lattice lie on the unit circle in the 
thermodynamic limit [|45j. Chen et al. [f|T] conjectured that when Q reaches a certain critical 



value Q C (L), all Fisher zeros for Lx L square lattices with self-dual boundary conditions are 
located at the unit circle \p\ = 1. In this paper we verify this conjecture and find that Q C {L) 
approaches infinity in the thermodynamic limit, and we study the thermal exponent y t of 
the square-lattice Potts antiferromagnet using the Fisher zeros near the antiferromagnetic 
critical point. 

In this paper we also discuss the partition function zeros in the complex Q plane (Potts 



zeros) of the Q-state Potts model. The Potts zeros at (3 J = — oo have been investigated 
extensively to understand the ground states of the antiferromagnetic Potts model and the 
chromatic polynomial in graph theory p3| , p6| , ^7| -[5T[1 . Recently the Potts zeros at finite 
temperatures have been studied for cyclic ladder graphs and e^ J < 1 [ |50|j . 

In the next section we describe two algorithms to evaluate the density of states, from 
which the exact partition function of the Q-state Potts model is obtained. The first algorithm 
(microcanonical transfer matrix) is applied to only integer Q but allows us to calculate the 
density of states for relatively larger lattices, while the second algorithm (random-cluster 
transfer matrix) gives the density of states for any value of Q. In Sec. Ill we discuss the 
Potts model at a = e^ J = 1 ± y/Q, its Potts zeros, and the related properties of the Fisher 
zeros. In the subsequent two sections we study the Potts zeros for the antiferromagnetic 
interval < a < 1 (Sec. IV) and for the ferromagnetic interval a > 1 (Sec. V). In Sec. VI 
we discuss the thermal exponent yt of the square lattice Q-state Potts antiferromagnet for 
2 < Q < 3 using the Fisher zeros. 



II. DENSITY OF STATES 

The Q-state Potts model for integer Q on a lattice G with N s sites and N b bonds is 
defined by the Hamiltonian 

where J is the coupling constant, (i, j) indicates a sum over nearest-neighbor pairs, 5 is the 
Kronecker delta, and crj = 1, 2, Q. The partition function of the model is 

Z Q =i: e~ m \ (2) 

where {a n } denotes a sum over Q Na possible spin configurations and (3 = (fc^T) -1 . If we 
define the density of states with energy < E < N b by 

fig(E)= E^-E^^i)). ( 3 ) 
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which takes on only integer values, then the partition function can be written as 

Z Q (a) = £ n Q (E)a E , (4) 

E=0 

where a = e^ J and states with E = (E = Nf,) correspond to the antiferromagnetic 
(ferromagnetic) ground states. From Eq. (4) it is clear that Zq(o) is simply a polynomial 
in a. We have calculated exact integer values for Qq =3 (E) of the three-state Potts model on 
finite L x L square lattices up to L = 12 using the microcanonical transfer matrix (/xTM) 

Here we describe briefly the /xTM |52| on an L x N square lattice with periodic boundary 
conditions in the horizontal direction (length L) and free boundaries in the vertical direction 
(length N). First, an array, w' 1 ', which is indexed by the energy E and variables <Tj, 1 < i < L 
for the first row of sites is initialized as 

L 

uW{E; a 2 , a L ) =5{E-J2 ^+i))- (5) 

i=l 

Now each spin in the row is traced over in turn, introducing a new spin variable from the 
next row, 

u(E; a[, a 2 , a L ) = u {1) (E - 8{a[, a{); a u a 2 , a L ). (6) 

This procedure is repeated until all the spins in the first row have been traced over, leaving 
a new function of the L spins in the second row. The horizontal bonds connecting the spins 
in the second row are then taken into account by shifting the energy, 

J 2 \E; a[, a> 2 , a' L ) = u(E - £ a[, a' 2 , a' L ). (7) 

i=l 

This procedure is then applied to each row in turn until the final (iVth) row is reached. The 
density of states is then given by 

W=EE-E " iN \E; a[, a' 2 , a' L ). (8) 

The permutation symmetry of the Q-state Potts model allows us to freeze the last spin 
<jl = 1 of each row. Now we need to consider only Q L ~ X possible spin configurations in each 
row instead of Q L configurations, and we save a great amount of memory and CPU time. 
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On the other hand, Fortuin and Kasteleyn |53| have shown that the partition function 
is also given by 

Z(a,Q)= £(a-l) fc ( G V (G,) , (9) 

G'CG 

where the summation is taken over all subgraphs G' C G, and b(G') and c(G') are, respec- 
tively, the number of occupied bonds and clusters in G'. In Eq. (9) Q need not be an 
integer and Eq. (9) defines the partition function of the (J-state Potts model for continuous 
Q. The random-cluster (or Fortuin-Kasteleyn) representation of the Potts model, Eq. (9), 
is also known as the Tutte dichromatic polynomial or the Whitney rank function in graph 
theory f50y5lfl . Introducing the density of states indexed by the number of occupied bonds 



< b < Nf, and the number of clusters 1 < c < N s , 

$(6, c) = £ 5(b - b{G'))5{c - c(G')), (10) 

G' 

which also takes on only integer values, the random-cluster representation of the Potts model 
can be written as 

6=0 c=l 

which is again a polynomial in a — 1 and Q. We have evaluated exact integer values for 
$(&, c) on finite LxL square lattices up to L = 8 for free, cylindrical, and self-dual boundary 
conditions using the random-cluster tranfer matrix. The self-dual lattices considered in this 
paper are periodic in the horizontal direction and there is another site above the LxL 
square lattice, which connects to L sites on the last (Lth) row (Figure 1). 

The algorithm (random-cluster transfer matrix) used to obtain the density of states 
$(6, c) is similar in spirit to that of Chen and Hu ||54|| . We consider anixiV square lattice 
with periodic boundary conditions in the horizontal direction (length L) and free boundaries 
in the vertical direction (length N). We define (p^ (b, c, {£}) as the density of states for the 
Lxm square lattice without the horizontal bonds in the mth row as a function of the number 
of occupied bonds b = 0, 1, 2L{m — 1), the number of clusters c = 1, 2, Lm, and the 



top labels {£} = {t\,t2, which tell whether each site in the mth row is connected to 

the other sites in the same row. 

The first step is to calculate <j>^ (6, c, {£}) using the Hoshen-Kopelman (HK) algorithm 



55| . The sites in the first row are labeled 1,2, ...,L from left to right and L + 1 to 2L in 
the second row. Cluster labels (i = 1, 2, 2L) are determined for each site and the top 
label tj (j = 1,2, ...,L) for the site j '• + L in the second row for each bond configuration. 
The top label tj is the smallest number of the set of indices j = 1, 2, L for the sites in the 
second row belonging to the same cluster which includes the site j + L. Because tj < j, the 
maximum number of sets of top labels {t} is L\. Counting the cases Sj = i gives the number 
of clusters c. 

Given (^^(b, c, {£}), 0^ m+1 ' ) (6, c, {t}) is calculated recursively by 
0i m+1) (6,c,{t}) = 7 E E 4> { r\b\d,{t!})5{b-b'-b g )5{c-d-^c)5{{t'}^{t}) (12) 

7=1 b',c',{t'} 

for m = 2, 3, N — 1, where 7 labels the 2 2L (= ^ max ) possible bond configurations in the 
newly added piece, g, consisting of the horizontal bonds in the mth row and the vertical 
bonds between the mth row and the (m + l)th row, and b g is the number of occupied bonds 
in g. The sites i in g are labeled from left to right by 1,2, ...,L in the mth row and by 
L + 1, L + 2, 2L in the (m + l)th row. We again use the HK algorithm to determine the 
cluster labels {si, S2, S2l} and the number of clusters c g in g, and the updated old top 
labels {?} = {^,£3, ■■■^'l} an d the new top labels {t} = {ti,t 2 , .-.^l} making a comparison 
between the cluster labels {si, s 2 , and the old top labels {£'} = {t'^t^ ■■■,t' L \. In Eq. 



'12) Ac is given by the Chen-Hu formula |p4 



Ac = c 9 -w-n' + n", (13) 

where n is the number of the cluster labels satisfying Sj = i for i — 1, 2, L, n/ the number 
of the old top labels satisfying t£ = z, and n" the number of the updated old top labels 
satisfying V i = i. 

Finally, the density of states $(6, c) is obtained by 
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Hb, c) = E E WW, WW -V- b g )5{c -d- Ac) 



(14) 



7=1 b',c',{t'} 



with jmax = 2 L and g made up of the horizontal bonds in the last (iVth) row. 

The random-cluster transfer matrix works very well, but for comparatively large lattices 
a considerable amount of memory is required to store (/>iP(b, c,{t}). At the expense of a 
slight increase in the complexity of the code it is possible to reduce the memory requirements 
substantially. First, the L\ sets of top labels include many unused sets, such as {...,U = 
i, tj = i,tk = j, ■■■} (i < j < k), which account for 56.7 % of all sets for L = 5 and 96.8 % for 
L — 10 and can be removed easily from 0^(6, c, {t}). Second, we should consider the fact 
that only some range of c is used for a fixed b. For example, in 4>f\b, c, {t}) only c = 1 to 
11 (5c = 11) are needed for b = 24. Here c = 1 results from the sparsest distributions of 24 
occupied bonds and c = 11 from the most compact distributions. 8c < 11 for all b ^ 24, and 
{8c) max = 11. We can calculate (5c) max easily for <f)^\b,c, {£}) and reduce a large amount 
of memory. Third, $(6, c) can be obtained directly from ^ m) (6', d , {t'}) (m < N - 1) 
with 'jmax = 2 L+2L ( iV ~ m ) using Eq. (14). This method decreases memory requirements but 
increases CPU time, while the former two methods reduce both the memory and CPU time 
requirements. In general, the random-cluster transfer matrix based on Eq. (12) is very fast, 
taking just 30 seconds on a PC with one PENTIUM 100 MHz CPU to obtain $(6, c) on the 
5x5 square lattice with free boundary conditions. 

The density of states £Iq(E) is related to the density of states $(6, c) by 



b=Ec=l v-v 

for integer Q. In Eq. (15) Q need not be an integer and Eq. (15) defines the density of 
states CIq(E) of the Q-state Potts model for noninteger Q. 




(15) 



III. THE CRITICAL POTTS MODEL 



At the ferromagnetic critical point, a c = 
Potts model becomes 



1 + -JQ, the partition function of the Q-state 
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/ b+2c 

Zcp = Y,®M[\[Q) > ( 16 ) 

b,c V ' 

which is a polynomial in \fQ. This defines what we refer to as the critical Potts model. 
Since b — N s — 1, c = 1 and 6 = 0, c = N s set the lowest and highest orders, respectively, in 
the polynomial, we can write Eq. (16) as 

Zcp = w Ns+1 E K r w r , (17) 

r=0 

where w = \[Q. The coefficients K r of the new polynomial Zqp satisfy 

N s -1 

J2 K r = 2^ (18) 
r=0 

and 

Ns-l 

E K r (-iy = o. (19) 

r=Q 

Table I shows the coefficients for the 8x8 square lattice with free boundary conditions. 

In addition to the ferromagnetic critical point a c = 1 + y/Q, the point d c = 1 — y/Q, 
which is sometimes referred to the unphysical critical point, also maps into itself under the 
dual transformation (a — 1) (a — 1) = Q [|l|. This leads us to consider the corresponding 
critical Potts partition function 

N s -1 

Z C p = w Ns+1 E K rv r , (20) 

r=0 

where w = —w. Evidently Z C p can be obtained from Z C p simply by continuing w to 
negative values. With this understanding we consider Zcp{w) for arbitrary complex values 
of w. Note that the map of the complex w plane on to the complex Q plane is now 2-to-l. 

Figure 2 shows the Potts zeros in the complex w plane of the critical Potts model on 
an 8 x 8 square lattice with self-dual boundary conditions. The zero at w = is N s + 1 
degenerate, and most of the remaining N s — 1 zeros lie in the half space Re(w) < 0. Several 
of these zeros lie on the negative real axis, and these will map on to the positive real Q axis 
as shown in Figure 3. Some of these zeros (Table II) lie at or close to the Beraha numbers 
47] 
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B n = 4cos 2 - (21) 
n 

with n = 2,3, 4, ... and < B n < 4. In the study of the phase diagram of the Potts model 
Saleur |H| assumed that the Potts model at the unphysical critical point, a c = 1 — \[Q, is 



singular when Q = B n , and our results verify this observation. Table II shows the Potts zeros 
of the critical Potts model on the L x L square lattice which lie at or close to the Beraha 
numbers for free (-/V& = 2L 2 — 2L), cylindrical (2L 2 — L), and self-dual (2L 2 ) boundary 
conditions. As the number of bonds, Nb, increases, the number of the Potts zeros at or close 
to the Beraha numbers B n increases for a fixed L, and as L increases the number of the 
Potts zeros at or close to B n increases for any specified type of boundary conditions. We 
expect that in the thermodynamic limit the Potts zeros on the positive real axis cover all 
the Beraha numbers B n (n = 2, 3, ...). 

For self-dual boundary conditions there exist unexpected Potts zeros on the positive real 
axis for Q > 4 (Table III). These zeros do not exist for non-dual boundary conditions, and 
the largest of these zeros, which we shall denote by Q ma x(L), has an interesting significance. 
Recently the partition function zeros in the complex temperature plane (Fisher zeros) have 
been studied extensively for the Potts model [2^,33-4^]. By numerical methods it has been 



shown [27,35,36,|0,|T,|3,53,|6| that for self-dual boundary conditions the Fisher zeros of the 
Q > 1 Potts models on a finite square lattice are located on the unit circle in the complex 
p plane for Re(p) > 0, where p = (a — l)/\/Q. Chen et al. [4"T| conjectured that when Q 



reaches a certain critical value Q C (L), all Fisher zeros are located on the unit circle \p\ = 1. 
However, the value of Q C (L) and how it scales with L were not addressed. We find that 
Q C {L) is identical to Q m ax{L) and that Q C {L) increases with L as shown in Table III. 

Figure 4 shows the Fisher zeros in the complex p plane of the Q-state Potts model on 
the 4x4 square lattice with self-dual boundary conditions. For Q = 75 the two zeros on 
the negative real axis lie off the unit circle, while for Q = 76 all the Fisher zeros lie on the 
unit circle. At Q = Q c (= 75.37... for L = 4) the two zeros lie on p = —1. In general, for 
the values of Q (both Q < 4 and Q > 4) that are determined from the Potts zeros on the 
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positive real axis, two Fisher zeros always lie at p = — 1. Q = 1 is exceptional in that all 
Fisher zeros of the one-state Potts model lie at p — — 1 [^TJ. Note that in Figure 4(b) the 
Fisher zeros are grouped and there exists a wide gap between two neighboring groups except 
for p = —1. Whenever all Fisher zeros lie on the unit circle, the number of groups of zeros 
is 2L X and the number of zeros for each group is L y , where L x and L y are the lattice sizes 
in the horizontal and vertical directions, respectively. 

By using the Bulirsch-Stoer (BST) algorithm ]5£]] we extrapolated l/Q c (L) for finite 
lattices to infinite size. The error estimates are twice the difference between the (n — 1, 1) 
and (n — 1, 2) approximants. For u — 1 (the parameter of the BST algorithm) we get 1/Q C = 
0.0007(8) and 1/Q C = 0.0001(7) for uj = 2. These results imply that in the thermodynamic 
limit all the Fisher zeros lie on the unit circle only in the limit Q — > oo ||42|| . Conversely, 
this observation implies that the locus of zeros in the thermodynamic limit for finite Q is 
an open question. 

IV. ANTIFERROMAGNETIC POTTS ZEROS 

For antiferromagnetic interaction, J < 0, the physical interval is0<a<l(0<T< oo). 
At zero temperature (a = 0) the partition function is 

Z = 5>(&,c)(-l) b Q c , (22) 

b,c 



which is also known as the chromatic polynomial in graph theory |5^j5^]. Figure 5 shows 



the zeros of the chromatic polynomial in the complex Q plane for the 8x8 square lattice for 
cylindrical |47] and self-dual boundary conditions. In Figure 5, except for the zeros at the 



Beraha numbers 0, 1 and 2 (Q = 2.0000000000007 for cylindrical boundary conditions), the 
Potts zeros are distributed along curves which cut the positive real axis between Q = 2 and 
3. The intersection of the locus of the Potts zeros with the real axis depends on the boundary 
condition: for L = 8 and cylindrical boundary conditions we have Q = 2.551073, while for 
self-dual boundary conditions we find a pair of zeros at Q = 2.636589 and 2.645969 which 
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are slightly larger than the fifth Beraha number B§ = 2.618034. For the 7x7 self-dual lattice 
these zeros lie at Q = 2.621577 and Q = 2.684634 (Figure 6). In addition for L = 7 there 
are isolated zeros on the real axis at the Beraha numbers B2 = 0, B3 = 1, and B4 = 2, and 
an additional zero appears at £> 6 = 3 (Figure 6). Q = 3 corresponds to the critical value Q c 
4"8| , |4Tj|] which separates the region (Q < 3) with antiferromagnetically ordered ground-states 



from the region (Q > 3) of disordered states at T = 0. Here we generalize this concept to 
finite temperatures and define Q c {a) to be the value of Q for a given value of a below which 
there is antiferromagnetic order. Because four colors are needed to color an L x L square 
lattice with self-dual boundary conditions such that no two nearest neighbors have the same 
color, there exists a trivial Potts zero at Q c = 3 when L = 3, 5, 7, ... . 

Figure 6 shows the Potts zeros of the dichromatic polynomial at several temperatures 
for the 7x7 lattice with self-dual boundary conditions. As a is increased the zeros move 
toward the origin and converge on the point Q = for a = 1 ||50|| . The antiferromagnetic 



critical point is given by a c (Q) = a/4 — Q — 1 fllUpTI , from which we have 

Q c (a) = (l-o)(a + 3). (23) 

Table IV shows the Potts zeros Q C (L) on (L = 4, 6, 8) or closest to (L = 3, 5, 7) the positive 
real axis for a = 0.5. From the BST extrapolation we obtained Q c = 1.78(18) (from 
L = 4,6,8) and Q c = 1.77(36) - 0.01(3)i (from L = 3,5,7) in agreement with Eq. (23). 
Figure 7 compares Eq. (23) (continuous curve) with the BST estimates from Q c (a, L) for 
L = 3, 5, 7 and self-dual boundary conditions for several values of a. 



V. FERROMAGNETIC POTTS ZEROS 

For ferromagnetic interaction, J > 0, the physical interval is a = [1, 00] (T = [00, 0]). 
Figure 8 shows the Potts zeros of the dichromatic polynomial on LxL lattices with cylindrical 
boundary conditions for a = 1 + \/2 = 2.414... and a = 1 + y/3 = 2.732... . For free and self- 
dual boundary conditions the distribution of the Potts zeros is similar to that for cylindrical 
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boundary conditions. Unlike the antiferromagnetic Potts zeros which are distributed mainly 
in the Ke(Q) > region (Figures 5 and 6), many ferromagnetic Potts zeros lie in the 
Ke(Q) < region. With the exception of the trivial zero at Q = the ferromagnetic Potts 
zeros are distributed along a single curve which moves away from the origin as a increases. 
There is no zero on the positive real axis, but the zero Qi(a, L) closest to the positive real 



axis approaches the real axis as L increases. As in the Yang-Lee theory |28|1 , we expect 
Qi(a, L) — > Q c (a) in the limit L — > oo. Table V shows the BST estimates from Qi(a, L) at 
a = 1 + a/2 and 1 + a/3 for different boundary conditions, suggesting that the locus of the 
Potts zeros cuts the positive real axis at Q c = 2 and 3, respectively, in the thermodynamic 
limit. From the ferromagnetic critical point, a c (Q) = 1 + y/Q, we obtain 

Q c (a) = (a - l) 2 , (24) 

which we have confirmed for a = 1 + a/2 and 1 + a/3 and other values of a > 1 (Figure 9). 
The behavior of the closest zero Q±(a, L) suggests a new scaling exponent y q defined as 

Q l {a,L)~Q c {a) + AL~ y ". (25) 

For finite lattices we define [FMBRlSlipj, 



Vq[L) ~ ln[(L + l)/L] • (26) 

The exponent y q is to the Potts zeros in the complex Q plane what the thermal exponent yt 
(or the magnetic exponent yn) is to the Fisher zeros in the complex temperature plane (the 
Yang-Lee zeros in the complex magnetic-field plane). Figure 10 shows the BST estimates 
from y q (L) for a = 2 (Q c = 1), 1 + a/2 (Q c = 2), 1 + a/3 (Q c = 3), and 3 (Q c = 4). 
The exponent y q increases as a (or Q c ) increases. Figure 10 compares our results for y q 
versus a c (Q) with the den Nijs formula @J2^] for the thermal exponent yt{a c {Q)) of the 
ferromagnetic Potts model. Clearly the general behaviors of y q and yt with a c (Q) are similar; 
these initial results are of insufficient precision to settle the question that y q = y t or not. 
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VI. FISHER ZEROS AND POTTS ANTIFERROMAGNETS 



For antiferromagnetic interaction J < the physical interval is < a = e f3J < 1 
(0 < T < oo), which corresponds to 

' < v = a -^r < o. (27) 



VQ~ VQ 

From the exact partition functions, Eqs. (4) and (11), we have evaluated Fisher zeros of 
the Potts model. Figure 11 shows the Fisher zeros in the complex p plane of the three-state 
Potts model on a 12 x 12 square lattice with free boundary conditions. The Fisher zeros in 
the complex p plane of the Q-state Potts model for several values of non-integer Q have been 
shown for the 8x8 square lattice with self-dual boundary conditions |27[] . Figure 12 shows the 
Fisher zeros in the complex p plane of the Q = 2.5 Potts model on an 8 x 8 square lattice with 
free boundary conditions. In Figures 11 and 12 there is a group of complex zeros approaching 
the antiferromagnetic critical point a c = \JA — Q — 1, equivalently, p c = (a c — and 
crossing the real axis at this critical point in the thermodynamic limit p7| . For an L x L 



square lattice a c (L) or p c (L) denotes the closest zero to the antiferromagnetic critical point 
or edge singularity. Based on the finite-size scaling law of the partition function zeros near 
the critical point |57|,58] we expect 



Im[a c (L)] ~ L~ Vt , (28) 

from which we can estimate the thermal exponent yt{L) for finite lattices as p7|j38|j39|j46|j52] 

ln{Im[q e (L + l)]/ImK(L)]} 
Vt{L) ~ ln[(L + I)/ L] • (29) 

Table VI shows the thermal exponents yt{L) of the Ising (Q = 2) antiferromagnet and the 

three-state Potts antiferromagnet for free boundary conditions. By using the BST algorithm 

we extrapolated our results for yt(L) to infinite size for 2 < Q < 3. Figure 13 shows the 

thermal exponent yt of the Potts antiferromagnet by the BST estimates with & — 1 (the 

parameter of the BST algorithm) for free boundary conditions. For the BST extrapolation 

of finite-size results of the Potts antiferromagnet we prefer free boundary condition to other 
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boundary conditions. The reason for this is that even though finite size effects are larger for 
free than cylindrical boundary conditions, the edge singularity approaches the critical point 
monotonically only if we consider a sequence of lattices with L even. For free boundary 
conditions this is not a problem and the increased effectiveness of the BST algorithm with 
longer sequences more than compensates the stronger finite-size effects [B3LB7[. In Figure 



13 there are two BST estimates for Q = 3. The upper estimate resulted from data for 
L = 3 ~ 8, while the lower one uses L = 3 ~ 12. In Figure 13 the continuous curve is the 
fit to the BST estimates with 

1 + Au + Bu 2 

Vt = C + Du ' (30) 

where 

2 ^VQ 

u = cos — — , (31) 

7T 2 

and A = -2.2821, B = -7.4390, C = 3.9818, and D = 7.4011. The variable u arises 
naturally in the expressions for the free energy f c (— |w) at the ferromagnetic and anti- 
ferromagnetic [[HJ critical temperatures, and in the critical exponents y t [0-0-0 and y h 



of the ferromagnetic Potts model. The form used in Eq. (30) has also been used to describe 
the critical exponent yh of the ferromagnetic Potts model Q. 

The BST estimates of the thermal exponent y t for Q < 3 are insensitive to the parameter 
of the BST algorithm, u>. However, as Q approaches 3 the BST results for the three-state 
Potts antiferromagnet are very sensitive to u. For example, we obtained y t = 0.50(8) for 
uj = \, y t = 0.59(2) for uj = 2, and y t = 0.60(2) for u = 3 using data for L = 3 ~ 12. 
The BST estimates of the thermal exponents of the Q-state Potts antiferromagnets for non- 



integer Q are also sensitive to u when Q ~ 3. Recently Ferreira and Sokal [20,24] have 



suggested the correlation length for the three-state Potts antiferromagnet has the form 

f ~ a - 1/yt {- lna) r (l + cia + c 2 a 2 + ...) (32) 
with yt = \ |T3J2nJ23], r w 1, and c\ ~ 15. For Q = 3 the sensitivity of the BST estimates 



of the thermal exponent to uj may result from this kind of logarithmic behavior. 
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Figure 14 shows the BST results extrapolated from Im[a c (L)] for L = 3 ~ 12 of the 
three-state Potts antiferromagnet with free boundary conditions as a function of uj along 
with the error estimates. When we use the BST algorithm to estimate a critical point, the 



best value of the free parameter to is the critical exponent yt We have obtained the 



desired result Im(a c ) = for uj = 0.5 which strongly suggests y t = 0.5. 

VII. CONCLUSION 

We have introduced the random-cluster transfer matrix to calculate exact integer values 
for the density of states $(6, c), from which the exact partition function Z(a,Q) can be 
obtained for any value of Q, even for complex Q. We have found a subset of the zeros of the 
partition function of the critical Potts model in the complex w = \[Q plane which lie close to 
or at the Beraha numbers on the negative real axis. The largest of these determines Q C (L), 
the value of Q above which the locus of zeros in the complex p plane lie on the unit circle. 
By studying the scaling behavior of Q C {L) with L we find that 1/Q C (L) — > as L —>■ oo, 
indicating that all the zeros do not lie strictly on the unit circle in the thermodynamic limit. 

We have studied the locus of zeros of the dichromatic polynomials in both the ferromag- 
netic and antiferromagnetic cases and find that the Yang-Lee mechanism is at work in the 
complex Q plane. We find Q^ F (a) = (1 — a) (a + 3) in agreement with Baxter |10,27|, and 
Qc M ( a ) = i a ~ I) 2 which is well known from duality arguments. Finally, we introduce a 
new finite-size scaling exponent, y q , which describes the approach of the edge singularity in 
the complex Q plane to the critical point as L — > oo. We find that y q varies with Q in much 
the same way as the thermal exponent y t of the ferromagnetic Potts model, but as yet we 
have not established a functional relation between y t and y q . 

We have also described the microcanonical transfer matrix to evaluate exact integer 
values for the density of states Qq(E) for the Q-state Potts model. From the densities 
of states $(6, c) and Qq(E) the partition functions Z(a,Q) and Zq(a) are obtained at any 
temperature a. Using the Fisher zeros of the exact partition functions we have estimated the 
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thermal exponents yt of the square-lattice Q-state Potts antiferro magnets for 2 < Q < 3. For 
Q < 3 the BST estimates are quite stable and y t is well approximated by a simple algebraic 
function of u — — ^cos -1 ^. However, as Q approaches 3, the BST estimates become 
sensitive to the choice of the scaling exponent u and to the data set used. Logarithmic 
or other corrections to scaling may be responsible for this behavior. For 3 < L < 8 and 
using the fit from data for Q < 3 we estimate yt{Q — 3) ~ 0.60(2), whereas if we include 
calculations for L up to 12 we find ?/t(3) ~ 0.50(8), in agreement with the leading scaling 
behavior suggested by Ferreira and Sokal [^0],|24| . We hope to resolve this issue by extending 
our exact calculations to larger lattices both exactly and by evaluating the density of states 
by microcanonical Monte Carlo sampling . 
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TABLES 

TABLE I. The coefficients K r of the partition function Zcp of the critical Potts model on the 
8x8 square lattice with free boundary conditions. 



K r 



K r 



126231322912498539682594816 1 

2 25524986518920425393717379072 3 

4 800610370286991686735405550336 5 

6 9347575153984981720573769774608 7 

8 54404758441262921869365590686720 9 

10 183329627865230663968273103188608 11 

12 391942582489345467968147273830784 13 

14 565568818070192070648821897874128 15 

16 576045479726330572980576680006144 17 

18 428419763789360447590812451240080 19 

20 238937966305748243499621822108592 21 

22 102094428513780610351844031072160 23 

24 34010305186209829834846747925664 25 

26 8960463658391600957849394069728 27 

28 1888880629020154547292686697440 29 

30 321508677911960109772525527808 31 

32 44483696316227122956909056000 33 

34 5024202380355112158475486704 35 

36 464007025651505425890675200 37 

38 35006515754308767635423136 39 

40 2149257909558929021370016 41 

42 106650313357232985654928 43 

44 4233470330438712180496 45 



2561398756299931321297272832 
166557700763955734137534296320 
3023834586769553668673015126432 
24326213916516119921387986971009 
106224421227588059984113069365972 
282506930412461406319413706064154 
492998772987796894034162031881014 
594803437106450324737629079389339 
515761419835859402146512922316166 
331188758886170694649818860535541 
161285868900631598864845612258887 
60729794216206721605782144017468 
17962439609348242109957007244868 
4227668735828771561070342983222 
800023985396669919928624375932 
122688547769932427716252035294 
15331317278052765348109117036 
1565743527537870861554921235 
1307342348007794922 1 1596986 
8903442105259073008726006 
491955405372613275069456 
21867081986237184782295 
772403311175092063841 
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46 132514803950430984480 

48 3208188678305076656 

50 58534057491001584 

52 776998275543312 

54 7144741728032 

56 42365906128 

58 144763280 

60 233296 

62 112 



47 21322374026497257618 

49 449814829279725547 

51 7036231117685951 

53 78304124284593 

55 584538167122 

57 2678567507 

59 6504139 

61 6265 

63 1 



TABLE II. The Potts zeros on the positive real Q axis of the critical Potts model which lie at 
or close to the Beraha numbers B n (n = 2, 3, ...). 



boundary condition 


free 


cylindrical 


self-dual 


self-dual 


system size 


8x8 


8x8 


5x5 


8x8 


B 2 = 














B 3 = 1 


1 


1 


1 


1 


B 4 = 2 


2.000000 


2.000000 


2.000000 


2.000000 


B 5 = 2.618034 


2.618034 


2.618034 


2.618055 


2.618034 


B 6 = 3 


3.000031 


3.000000 


2.992072 


3.000000 


B 7 = 3.246980 


3.226656 


3.246976 




3.246980 


B 8 = 3.414214 




3.415672 


3.412158 


3.414685 


B 9 = 3.532089 




3.521330 




3.524855 


Bio = 3.618034 








3.618701 


B ie = 3.847759 








3.839893 


B 30 = 3.956295 






3.957208 




B 64 = 3.990369 








3.990438 
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TABLE III. The Potts zeros on the positive real axis for Q > 4 for the L x L square lattice 
with self-dual boundary conditions. 

L=4 5 6 7 8 

75.373518 185.886317 395.130118 754.036414 1324.684018 

7.566911 21.911010 40.294754 66.309209 

6.401881 15.678097 

5.326082 



TABLE IV. The Potts zeros on or closest to the positive real axis for the L x L square lattice 
with self-dual boundary conditions at a = 0.5. 

L Q C (L) L Q C (L) 



3 1.279400 + 0.161071* 4 1.441800 

5 1.499871 +0.0695198* 6 1.574011 

7 1.583953 + 0.0407605* 8 1.632666 





TABLE V. The BST estimate 


s from Q\(a,L) for different 


boundary conditions. 


a 


free 


cylindrical 


self-dual 


1 + V2 


1.90(10) +0.09(24)i 


1.94(7) +0.15(14)* 


1.95(10) +0.13(11)* 


i + Vs 


2.84(9) - 0.13(31)* 


2.88(2) +0.00(15)* 


2.89(8) - 0.03(10)*' 
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TABLE VI. The thermal exponents yt(L) of the Q-state Potts antiferromagnets for Q = 2 and 



Q = 


3 with free boundary conditions. The last row is the BST extrapolation with 


lo = 1 to infinite 


size. 






L 


Vt(L) (Q = 2) 


y t (L) (Q = 3) 


3 


0.859670530424 


0.672417300113 


4 


0.882900616441 


0.840771366429 


5 


0.895500892567 


0.750192805568 


6 


0.904846051999 


0.714132507277 


7 


0.912493138251 


0.694522575800 


8 


0.918981910221 


0.681414203729 


9 


0.924586147759 


0.671514256321 


10 


0.929481322004 


0.663473505003 


11 


0.933794047470 


0.656641075731 


oo 


1.000005(9) 


0.50(8) 
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FIGURES 




FIG. 1. 5 x 5 square lattice with self-dual boundary conditions. 
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FIG. 2. Potts zeros in the complex w (= Q 2 ) plane of the partition function Zqp for the 8 x 
square lattice with self-dual boundary conditions. 
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FIG. 3. Potts zeros in the complex Q plane of the critical Potts model for the 8x8 square 
lattice with free boundary conditions. 
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FIG. 4. Fisher zeros in the complex p plane of the Q-state Potts model on the 4x4 square 
lattice with self-dual boundary conditions for (a) Q = 75 and (b) Q = 76. 
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FIG. 5. Potts zeros in the complex Q plane of the chromatic polynomial on the 8x8 square 
lattice for cylindrical [15] and self-dual boundary conditions. 
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FIG. 6. Potts zeros of the dichromatic polynomial for < a < 1 on the 7x7 square lattice 
with self-dual boundary conditions. 
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FIG. 7. BST extrapolation of Q^ F as a function of a for self-dual boundary conditions. The 
continuous curve is given by Q c = (1 — a) (a + 3). 
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FIG. 8. Potts zeros in the complex Q plane of the dichromatic polynomial on the L x L square 
lattices with cylindrical boundary conditions for a > 1. 
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FIG. 9. BST extrapolation of Qf M as a function of a for self-dual boundary conditions. The 
continuous curve is given by Q c = (a — l) 2 . 
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FIG. 10. The exponent y q as a function of a for free, cylindrical, and self-dual boundary 
conditions. The slight horizontal offset for data for cylindrical and self-dual boundary conditions 
is for clarity only. The long-dashed curve is the thermal exponent yt by the den Nijs formula. 



37 



1 .2 



0.4 



CL 



■0.4 



-1 .2 



AAA A A 



AAA 
J 

A A 



A A 



A 
A 



A A A A 



A A A 
AA 

AA , t 



X- 



.A 



a 12x1 2 lattice 
x critical point 



A a 



AAA 



A A 



A A 
A A A 



AAA A ^ 



4 a Aa A J 
^ AAAa Vaaaaaaaaaaaaaa^ aAA 



-1 .5 



-0.5 



0.5 



Re(p) 



FIG. 11. Fisher zeros in the complex p plane of the three-state Potts model on 12 x 12 square 
lattice with free boundary conditions. The dotted line is the antiferromagnetic interval. 
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FIG. 12. Fisher zeros in the complex p plane of the Q = 2.5 Potts model on 8 x 8 square lattice 
with free boundary conditions. The dotted line shows the antiferromagnetic interval. 
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FIG. 13. The thermal exponents yt of the Q-state Potts antiferromagnets by the BST estimates 
(filled circles) from data for L = 3 ~ 8 and free boundary conditions. For Q = 3 the BST estimate 
(filled triangle) from data for L = 3 ~ 12 is added and has the slight horizontal offset for clarity 
only. 
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FIG. 14. BST extrapolation of the imaginary part of the critical point, Im(a c ), for 
three-state Potts antiferromagnet as a function of the parameter oj. 
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